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23.1. Introduction 

Tile energy content in radiation from beyond our Galaxy is 
dominated by the Cosmic Microwave Background (CMB), discovered 
in 1965 [1]. The spectrum of the CMB is well described by a 
blackbody function with T = 2.725 K. This spectral form is one of the 
main pillars of the hot Big Bang model for the early Universe. The 
lack of any observed deviations from a blackbody spectrum constrains 
physical processes over cosmic history at redshifts z < 10 (see earlier 
versions of this review). All viable cosmological models predict a very 
nearly Planckian spectrum inside the current observational limits. 

Another observable quantity inherent in the CMB is the variation 
in temperature (or intensity) from one part of the microwave sky to 
another [2]. Since the first detection of these anisotropies by the 
COBE satellite [3] , there has been intense activity to map the sky at 
increasing levels of sensitivity and angular resolution by ground-based 
and balloon-borne measurements. These were joined in 2003 by the 
first results from NASA's Wilkinson Microwave Anisotropy Probe 
(WMAP) [4], which were improved upon by analysis of the 3 
year and 5 year WMAP data [5,6]. Together these observations 
have led to a stunning confirmation of the 'Standard Model of 
Cosmology.' In combination with other astrophysical data, the CMB 
anisotropy measmements place quite precise constraints on a immber 
of cosmological parameters, and have launched us into an era of 
precision cosmology. This is expected to continue with the improved 
capabilities of the Planck satellite. 

23.2. Description of CMB Anisotropies 

Observations show that the CMB contains anisotropies at the 
10^° level, over a wide range of angular scales. These anisotropies 
are usually expressed by using a spherical harmonic expansion of the 
CMB sky: 

T(0,4>)=Y,ai,„Yi„,{B,4,). 

The vast majority of the cosmological information is contained in 
the temperature 2-point function, i.e., the variance as a function 
only of angular separation, since we notice no preferred direction. 

Equivalently, the power per unit In^ is iJ2m l^-tn 
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23.2.1. The Monopole : 

The CMB has a mean temperature of T^ = 2.725 ± 0.001 K (Icr) [7], 
which can be considered as the monopole component of CMB maps, 
ago- Since all mapping experiments involve difference measurements, 
they are insensitive to this average level. Monopole measurements 
can only be made with absolute temperature devices, such as the 
FIRAS instrument on the COBE satellite [7]. Such measurements 
of the spectrum are consistent with a blackbody distribution over 
more than three decades in frequency (with some recent evidence 
for deviation at low frequencies [8]) . A blackbody of the measured 
temperature corresponds to n-y = (2f (3)/7r ) T.^ ~ 411 cm^ and 
p^ = (7r2/15)T* ~ 4.64 x 10^34 g ^.j^-S ^ o.260eVcm-3, 

23.2.2. The Dipole : 

The largest anisotropy is in the t. = I (dipole) first spherical 
harmonic, with amplitude 3.355 ± 0.008 niK [6]. The dipole is 
interpreted to be the result of the Doppler shift caused by the solar 
system motion relative to the nearly isotropic blackbody field, as 
confirmed by measurements of the radial velocities of local galaxies [9] . 
The motion of an observer with velocity /3 = v/c relative to an 
isotropic Planckian radiation field of temperature Tq produces a 
Doppler-shifted temperature pattern 

T(e) = To (1-/32)'^'/ (1-/3 cose) 

~ro (l + l3cose+ (b'^ 1 2\ COS2B + O (l3'^y\ . 

At every point in the sky, one observes a blackbody spectrum, with 
temperature T(d). The spectrum of the dipole is the differential of a 
blackbody spectrum, as confirmed by Rcf. 10. 



The implied velocity for the solar system barycenter is v = 
369.0 ± 0.9 km s^, assuming a value Tq = T-y, towards {i,b) = 
(263.99° ± 0.14°, 48.26° ± 0.03°) [6,11]. Such a solar system motion 
implies a velocity for the Galaxy and the Local Group of galaxies 
relative to the CMB. The derived value is i'lq = 627±22kms^ 
towards {I, b) = (276° ± 3°, 30° ± 3°), where most of the error comes 
from uncertainty in the velocity of the solar system relative to the 
Local Group. 

The dipole is a frame-dependent quantity, and one can thus 
determine the 'absolute rest frame' as that in which the CMB dipole 
would be zero. Our velocity relative to the Local Group, as well as the 
velocity of the Earth around the Sun, and any velocity of the receiver 
relative to the Earth, is normally removed for the purposes of CMB 
anisotropy study. 

23.2.3. Higher- Order Multipoles : 

The variations in the CMB temperature maps at higher multipoles 
{I > 2) are interpreted as being mostly the result of perturbations 
in the density of the early Universe, manifesting themselves at the 
epoch of the last scattering of the CMB photons. In the hot Big Bang 
picture, the expansion of the Universe cools the plasma so that by a 
redshift z ~ 1100 (with little dependence on the details of the model), 
the hydrogen and helium nuclei can bind electrons into neutral atoms, 
a process usually referred to as recombination [12]. Before this epoch, 
the CMB photons are tightly coupled to the baryons, while afterwards 
they can freely stream towards us. 

Theoretical models generally predict that the a(^ modes are 
Gaussian random fields to high precision, e.g., standard slow-roll 
inflation's non-Gaussian contribution is expected to be one or two 
orders of magnitude below current observational limits [13] . Although 
non-Gaussianity of various forms is possible in early Universe models, 
tests show that Gaussianity is an extremely good simplifying 
approximation [14], with only some relatively weak indications 
of non-Gaussianity or statistical anisotropy at large scales. Such 
signatures found in existing WMAP data are generally considered to 
be subtle foreground or instrumental artefacts [15,16]. 

A statistically isotropic sky means that all ms are equivalent, i.e., 
there is no preferred axis. Together with the assumption of Gaussian 
statistics, the variance of the temperature field (or equivalently the 
power spectrum in £) then fully characterizes the anisotropies. The 
power summed over all ms at each I is {2i + 1)C(/{4tt), where 
Ci = (jafml )• Thus averages of ag^s over m can be used as 
estimators of the CfS to constrain their expectation values, which 
are the quantities predicted by a theoretical model. For an idealized 
full-sky observation, the variance of each measured Ci {i.e., the 
variance of the variance) is [2/(2^ + 1)]C|. This sampling uncertainty 
(known as 'cosmic variance') comes about because each Cg is x 
distributed with {2( + 1) degrees of freedom for our observable volume 
of the Universe. For fractional sky coverage, /sky, this variance is 
increased by l/Zsky Eind the modes become partially correlated. 

It is important to understand that theories predict the expectation 
value of the power spectrum, whereas our sky is a single realization. 
Hence the cosmic variance is an unavoidable source of uncertainty 
when constraining models; it dominates the scatter at lower £s, while 
the effects of instrumental noise and resolution dominate at higher ^s 
[17]. 

23.2.4. Angular Resolution and Binning : 

There is no one-to-one conversion between nmltipole I and the 
angle subtended by a particular spatial scale projected onto the sky. 
However, a single spherical harmonic Yi^ corresponds to angular 
variations of S ~ tt/1. CMB maps contain anisotropy information from 
the size of the map (or in practice some fraction of that size) down 
to the beam-size of the instrument, a. One can think of the effect of 
a Gaussian beam as rolling off' the power spectrum with the function 

For less than full sky coverage, the i modes become correlated. 
Hence, experimental results are usually quoted as a series of 'band 
powers', defined as estimators of i{£ + l)Ci/2w over different ranges 
of £. Because of the strong foreground signals in the Galactic Plane, 
even 'all-sky' surveys, such as COBE and WMAP involve a cut sky. 
The amount of binning required to obtain uncorrelated estimates of 
power also depends on the map size. 



2 23. Cosmic microwave background 



23.3. Cosmological Parameters 

The current 'Standard Model' of cosmology contains around 10 
free parameters (see The Cosmological Parameters — Sec. 21 of this 
Review). The basic framework is the Friedmann- Robertson- Walker 
(FRW) metric (i.e., a universe that is approximately homogeneous 
and isotropic on large scales), with density perturbations laid down 
at early times and evolving into today's structures (see Big-Bang 
cosmology — Sec. 19 of this Review). The most general possible set of 
density perturbations is a linear combination of an adiabatic density 
perturbation and some isocurvature perturbations. Adiabatic means 
that there is no change to the entropy per particle for each species, 
i.e., &pI p for matter is {3/4)5p/p for radiation. Isocurvature means 
that the set of individual density perturbations adds to zero, for 
example, matter perturbations compensate radiation perturbations 
so that the total energy density remains unperturbed, i.e., 5p for 
matter is —Sp for radiation. These different modes give rise to 
distinct (temporal) phases during growth, with those of the adiabatic 
scenario being strongly preferred by the data. Models that generate 
mainly isocurvature type perturbations (such as most topological 
defect scenarios) are no longer considered to be viable. However, an 
admixture of the adiabatic mode with up to about 10% isocurvature 
contribution is still allowed [18]. 

Within the adiabatic family of models, there is, in principle, a 
free function describing how the comoving curvature perturbations, 
7^(x, t) , vary with length scale. The great virtue of TZ is that it 
is constant for a purely adiabatic perturbation. There are physical 
reasons to anticipate that the variance of these perturbations will 
be described well by a power-law in scale, i.e., in Fourier space 
(|7^||) oc fe"^ , where k is wavenumber and n is the usual definition of 
spectral index. So-called 'scale-invariant' initial conditions (meaning 
gravitational potential fiuctuations which are independent of k) 
correspond to n = 1. In inflationary models [19], perturbations are 
generated by quantum fluctuations, which are set by the energy scale 
of inflation, together with the slope and higher derivatives of the 
inflationary potential. One generally expects that the Taylor series 
expansion of In TZk (In fc) has terms of steadily decreasing size. For 
the simplest models, there are thus 2 parameters describing the 
initial conditions for density perturbations: the amplitude and slope 
of the power spectrum. These can be explicitly defined, for example, 
through: 



= (fcV2^') 
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with A = A^(fco) and fcg = 0.002 Mpc^, say. There arc many 
other equally valid definitions of the amplitude parameter (see also 
Sec. 19 and Sec. 21 of this Review), and we caution that the 
relationships between some of them can be cosmology-dependent. In 
'slow roll' inflationary models, this normalization is proportional to 
the combination V / {V') , for the inflationary potential V{tj>). The 
slope n also involves V" , and so the combination of A and n can, in 
principle, constrain potentials. 

Inflation generates tensor (gravitational wave) modes, as well as 
scalar (density perturbation) modes. This fact introduces another 
parameter, measuring the amplitude of a possible tensor component, or 
equivalently the ratio of the tensor to scalar contributions. The tensor 
amplitude is A'y oc V , and thus one expects a larger gravitational wave 
contribution in models where inflation happens at higher energies. 
The tensor power spectrum also has a slope, often denoted nx, but 
since this seems unlikely to be measured in the near future, it is 
sufficient for now to focus only on the amplitude of the gravitational 
wave component. It is most common to define the tensor contribution 
through r, the ratio of tensor to scalar perturbation spectra at some 
small value of k (although sometimes it is defined in terms of the 
ratio of contributions at ^ = 2). Different inflationary potentials will 
lead to different predictions, e.g., for \(j) inflation with 50 e- folds, 
r = 0.32, and for m <p inflation r ~ 0.15, while other models can 
have arbitrarily small values of r. In any case, whatever the specific 
definition, and whether they come from inflation or something else, 
the 'initial conditions' give rise to a minimum of 3 parameters: A, n, 
and r. 

The background cosmology requires an expansion parameter 
(the Hubble Constant, Hq, often represented through Hq = 
100/ikms^ Mpc^ ) and several parameters to describe the matter 
and energy content of the Universe. These are usually given in terms 
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Figure 23.1: The theoretical CMB anisotropy power spectrum, 
using a standard ACDM model from CMBFAST. The z-axis is 
logarithmic here. The regions, each covering roughly a decade in 
(., are labeled as in the text: the ISW rise; Sachs- Wolfe plateau; 
acoustic peaks; and damping tail. Also shown is the shape of 
the tensor (gravitational wave) contribution, with an arbitrary 
normalization. 

of the critical density, i.e., for species 'x'. Ox = Px/Pcriti where 
Pcrit = 3Hq/8ttG. Since physical densities px oc Qy^h = Wx are what 
govern the physics of the CMB anisotropics, it is these uis that are best 
constrained by CMB data. In particular CMB observations constrain 
Oj,/i for baryons and Hm/i for baryons plus Cold Dark Matter. 

The contribution of a cosmological constant A (or other form of 
Dark Energy) is usually included via a parameter which quantifies 
the curvature, Ok = 1 — fitoti where Otot = f^m + f^A- The radiation 
content, while in principle a free parameter, is precisely enough 
determined by the measurement of T-y, and makes a < 10^ 
contribution to fitot today. 

The main effect of astrophysical processes on the C'lS comes 
through reionization. The Universe became reionized at some redshift 
2;, long after recombination, affecting the CMB through the integrated 
Thomson scattering optical depth: 



"^ 



where irx is the Thomson cross-section, ne{z) is the number density 
of free electrons (which depends on astrophysics), and dt/dz is fixed 
by the background cosmology. In principle, r can be determined from 
the small-scale matter power spectrum, together with the physics 
of structure formation and feedback processes. However, this is a 
sufficiently intricate calculation that r needs to be considered as a free 
parameter. 

Thus, we have 8 basic cosmological parameters: A, n, r, h, Q.^h , 
Q.xnh , Ototi and r. One can add additional parameters to this list, 
particularly when using the CMB in combination with other data sets. 
The next most relevant ones might be: Cl^h , the massive neutrino 
contribution; w (= p/p), the equation of state parameter for the 
Dark Energy; and dn/d\nk, measuring deviations from a constant 
spectral index. To these 1 1 one could of course add further parameters 
describing additional physics, such as details of the reionization 
process, features in the initial power spectrum, a sub-dominant 
contribution of isocurvature modes, etc. 

As well as these underlying parameters, there are other quantities 
that can be obtained from them. Such derived parameters include the 
actual Os of the various components (e.g., fim), the variance of density 
perturbations at particular scales [e.g., ag,), the age of the Universe 
today (to) I the age of the Universe at recombination, reionization, etc. 
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23.4. Physics of Anisotropies 

The cosmological parameters affect the anisotropies through the 
well understood physics of the evolution of linear perturbations within 
a background FRW cosmology. There are very effective, fast, and 
publicly- available software codes for computing the CMB anisotropy, 
polarization, and matter power spectra, e.g., CMBFAST [20] and 
CAMB [21]. These have been tested over a wide range of cosmological 
parameters and are considered to be accurate to better than the 1% 
level [22]. 

A description of the physics underlying the C^s can be separated 
into 3 main regions, as shown in Fig. 23.1. 

23.4.1. The ISW rise, e ^ 10, and Sachs-Wolfe plateau, 
10 ^ £ ^ 100 : 

The horizon scale (or more precisely, the angle subtended by the 
Hubble radius) at last scattering corresponds to ^ ~ 100. Anisotropies 
at larger scales have not evolved significantly, and hence directly 
refiect the 'initial conditions'. 5T/T = -(1/5)7?.(xlss) ^ (l/3)5(/i/c2, 
here (50 is the perturbation to the gravitational potential, evaluated on 
the last scattering surface (LSS). This is a result of the combination 
of gravitational redshift and intrinsic temperature fluctuations and is 
usually referred to as the 'Sachs- Wolfe' effect [23]. 

Assuming that a nearly scale-invariant spectrum of curvature 
and corresponding density perturbations was laid down at early 
times {i.e., n ~ 1, meaning equal power per decade in k), then 
({( + l)Ci ~ constant at low £s. This effect is hard to see unless the 
multipole axis is plotted logarithmically (as in Fig. 23.1, but not 
Fig. 23.2). 

Time variation of the potentials (i.e., time-dependent metric 
perturbations) leads to an upturn in the C(S in the lowest several 
multipoles; any deviation from a total equation of state w = has 
such an effect. So the dominance of the Dark Energy at low redshift 
makes the lowest ^s rise above the plateau. This is sometimes called 
the 'integrated Sachs-Wolfe effect' (or ISW rise), since it comes from 
the line integral of <f>; it has been confirmed through correlations 
between the large-angle anisotropies and large-scale structure [24]. 
Specific models can also give additional contributions at low I 
{e.g., perturbations in the Dark Energy component itself [25]), but 
typically these are buried in the cosmic variance. 

In principle, the mechanism that produces primordial perturbations 
could generate scalar, vector, and tensor modes. However, the vector 
(vorticity) modes decay with the expansion of the Universe. The 
tensors (transverse trace-free perturbations to the metric) generate 
temperature anisotropies through the integrated effect of the locally 
anisotropic expansion of space. Since the tensor modes also redshift 
away after they enter the horizon, they contribute only to angular 
scales above about 1° (see Fig. 23.1). Hence some fraction of the low 
£ signal could be due to a gravitational wave contribution, although 
small amounts of tensors are essentially impossible to discriminate 
from other effects that might raise the level of the plateau. However, 
the tensors can be distinguished using polarization information (see 
Sec. 23.6). 

23.4.2. The acoustic peaks, 100 ^ £ ^ 1000 : 

On sub-degree scales, the rich structure in the anisotropy spectrum 
is the consequence of gravity-driven acoustic oscillations occurring 
before the atoms in the Universe became neutral. Perturbations 
inside the horizon at last scattering have been able to evolve causally 
and produce anisotropy at the last scattering epoch, which reflects 
this evolution. The frozen-in phases of these sound waves imprint 
a dependence on the cosmological parameters, which gives CMB 
anisotropies their great constraining power. 

The underlying physics can be understood as follows. Before the 
Universe became neutral, the proton-electron plasma was tightly 
coupled to the photons, and these components behaved as a single 
'photon-baryon fluid.' Perturbations in the gravitational potential, 
dominated by the Dark Matter component, were steadily evolving. 
They drove oscillations in the photon-baryon fluid, with photon 
pressure providing most of the restoring force and baryons giving some 
additional inertia. The perturbations were quite small in amplitude, 
0(10^ ), and so evolved linearly. That means each Fourier mode 
developed independently, and hence can be described by a driven 
harmonic oscillator, with frequency determined by the sound speed in 



the fluid. Thus the fluid density underwent oscillations, giving time 
variations in temperature. These combine with a velocity effect which 
is tt/2 out of phase and has its amplitude reduced by the sound speed. 

After the Universe recombined, the radiation decoupled from the 
baryons and could travel freely towards us. At that point, the phases 
of the oscillations were frozen-in, and became projected on the sky 
as a harmonic series of peaks. The main peak is the mode that went 
through 1/4 of a period, reaching maximal compression. The even 
peaks are maximal anrfer-densities, which are generally of smaller 
amplitude because the rebound has to fight against the baryon inertia. 
The troughs, which do not extend to zero power, are partially filled by 
the Doppler effect because they are at the velocity maxima. 

The physical length scale associated with the peaks is the sound 
horizon at last scattering, which can be straightforwardly calculated. 
This length is projected onto the sky, leading to an angular scale that 
depends on the geometry of space, as well as the distance to last 
scattering. Hence the angular position of the peaks is a sensitive probe 
of the spatial curvature of the Universe {i.e., Htot), with the peaks 
lying at higher I in open universes and lower £ in closed geometry. 

One additional effect arises from reionization at redshift 2;. A 
fraction of photons (r) will be isotropically scattered a.t z < z\, 
partially erasing the anisotropies at angular scales smaller than those 
subtended by the Hubble radius at Zj. This corresponds typically to ^s 
above about a few 10s, depending on the specific reionization model. 
The acoustic peaks are therefore reduced by a factor e^ relative to 
the plateau. 

These peaks were a clear theoretical prediction going back to about 
1970 [26]. One can think of them as a snapshot of stochastic standing 
waves. Since the physics governing them is simple and their structure 
rich, then one can see how they encode extractable information about 
the cosmological parameters. Their empirical existence started to 
become clear around 1994 [27], and the emergence, over the following 
decade, of a coherent series of acoustic peaks and troughs is a triumph 
of modern cosmology. This picture has received further confirmation 
with the detection in the power spectrum of galaxies (at redshifts 
close to zero) of the imprint of these same acoustic oscillations in the 
baryon component [28,29,30]. 

23.4.3. The damping tail, t ^ 1000 : 

The recombination process is not instantaneous, giving a thickness 
to the last scattering surface. This leads to a damping of the 
anisotropies at the highest ^s, corresponding to scales smaller 
than that subtended by this thickness. One can also think of the 
photon-baryon fiuid as having imperfect coupling, so that there is 
diffusion between the two components, and hence the amplitudes of 
the oscillations decrease with time. These effects lead to a damping 
of the Ci'A, sometimes called Silk damping [31], which cuts off the 
anisotropies at multipoles above about 2000. 

An extra effect at high ^s comes from gravitational lensing, caused 
mainly by non-linear structures at low redshift. The C^s are convolved 
with a smoothing function in a calculable way, partially fiattening 
the peaks, generating a power-law tail at the highest multipoles, and 
complicating the polarization signal [32] . The effects of lensing on 
the CMB have recently been detected by correlating temperature 
gradients and small-scale filtered anisotropies from WTsJlAV with 
lensing potentials traced using galaxies [33]. This is an example 
of a 'secondary effect,' i.e., the processing of anisotropies due to 
relatively nearby structures (see Sec. 23.7.2). Galaxies and clusters 
of galaxies give several such effects; all are expected to be of low 
amplitude and typically affect only the highest &, but they carry 
additional cosmological information and will be increasingly important 
as experiments push to higher sensitivity and angular resolution. 

23.5. Current Anisotropy Data 

There has been a steady improvement in the quality of CMB 
data that has led to the development of the present-day cosmological 
model. Probably the most robust constraints currently available come 
from the combination of the WMAP five year data [34] with smaller 
scale results from the ACBAR [35] and QUAD [36] experiments 
(together with constraints from other cosmological data-sets). We 
plot power spectrum estimates from these experiments, as well 
as BOOMERANG [38] and CBI [39] in Fig. 23.2. Other recent 
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experiments also give powerful constraints, which are quite consistent 
with what we describe below. There have been some comparisons 
among data-sets [37,41], which indicate very good agreement, both 
in maps and in derived power spectra (up to systematic uncertainties 
in the overall calibration for some experiments). This makes it clear 
that systematic effects are largely under control. However, a fully 
self-consistent joint analysis of all the current data sets has not 
been attempted, one of the reasons being that it requires a careful 
treatment of the overlapping sky coverage. 
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Figure 23.2: Band-power estimates from the WMAP, 
BOOMERANG, QUAD, CBI, and ACBAR experiments. 
Some of the low-i and high-f band-powers which have large 
error bars have been omitted. Note also that the widths of the 
£-bands varies between experiments and have not been plotted. 
This figure represent only a selection of available experimental 
results, with some other data-sets being of similar quality. The 
multipole axis here is linear, so the Sachs- Wolfe plateau is hard 
to see. However, the acoustic peaks and damping region are very 
clearly observed, with no need for a theoretical curve to guide 
the eye; the curve plotted is a best-fit model from WMAP 5-year 
plus other CMB data. 

The band-powers shown in Fig. 23.2 are in very good agreement 
with a 'ACDM' type model, as described earlier, with several of the 
peaks and troughs quite apparent. For details of how these estimates 
were arrived at, the strength of any correlations between band-powers 
and other information required to properly interpret them, the original 
papers should be consulted. 

23.6. CMB Polarization 

Since Thomson scattering of an anisotropic radiation field also 
generates linear polarization, the CMB is predicted to be polarized 
at the roughly 5% level of the temperature anisotropics [42]. 
Polarization is a spin-2 field on the sky, and the algebra of the modes 
in ^-space is strongly analogous to spin-orbit coupling in quantum 
mechanics [43]. The linear polarization pattern can be decomposed 
in a number of ways, with two quantities required for each pixel in 
a map, often given as the Q and U Stokes parameters. However, 
the most intuitive and physical decomposition is a geometrical one, 
splitting the polarization pattern into a part that comes from a 
divergence (often referred to as the 'E-mode') and a part with a curl 
(called the 'B-mode') [44]. More explicitly, the modes are defined in 
terms of second derivatives of the polarization amplitude, with the 
Hessian for the E-modes having principle axes in the same sense as 
the polarization, while the B-mode pattern can be thought of simply 
as a 45° rotation of the E-mode pattern. Globally one sees that the 
E-modes have (—1) parity (like the spherical harmonics), while the 
B-modes have (— 1) + parity. 

The existence of this linear polarization allows for 6 different 
cross power spectra to be determined from data that measure the 
full temperature and polarization anisotropy information. Parity 



considerations make 2 of these zero, and we are left with 4 potential 
observables: Cj , Cj , Cf , and C^ . Because scalar perturbations 
have no handedness, the B-mode power spectrum can only be sourced 
by vectors or tensors. Moreover, since inflationary scalar perturbations 
give only E-modes, while tensors generate roughly equal amounts of 
E- and B-modes, then the determination of a non-zero B-mode signal 
is a way to measure the gravitational wave contribution (and thus 
potentially derive the energy scale of inflation), even if it is rather 
weak. However, one must first eliminate the foreground contributions 
and other systematic effects down to very low levels. 

The oscillating photon-baryon fluid also results in a series of acoustic 
peaks in the polarization Qs. The main 'EE' power spectrum has 
peaks that are out of phase with those in the 'TT' spectrum, because 
the polarization anisotropics are sourced by the fluid velocity. The 
'TE' part of the polarization and temperature patterns comes from 
correlations between density and velocity perturbations on the last 
scattering surface, which can be both positive and negative, and is of 
larger amplitude than the EE signal. There is no polarization 'Sachs- 
Wolfe' effect, and hence no large-angle plateau. However, scattering 
during a recent period of reionization can create a polarization 'bump' 
at large angular scales. 

Because the polarization anisotropics have only a fraction of the 
amplitude of the temperature anisotropics, they took longer to detect. 
The first measurement of a polarization signal came in 2002 from 
the DASI experiment [45], which provided a convincing detection, 
confirming the general paradigm, but of low enough significance that 
it lent little constraint to models. As well as the E-mode signal, DASI 
also made a statistical detection of the TE correlation. 

In 2003, the WMAP experiment demonstrated that is was able 
to measure the TE cross-correlation power spectrum with high 
precision [46], and this was improved upon in the 3- and 5-year 
results, which also included EE measurements [34,47]. The TE signal 
has been mapped out quite accurately now through a combination 
of data from WMAP, together with the BOOMERANG [54], 
QUAD [37], CBI [51], DASI [52] and BICEP [49,50] experiments, 
which are shown in Fig. 23.3. The anti-correlation at ^ ~ 150 and the 
peak at ^ ~ 300 are now quite distinct. The measured shape of the 
cross-correlation power spectrum provides supporting evidence for the 
adiabatic nature of the perturbations, as well as directly constraining 
the thickness of the last scattering surface. Since the polarization 
anisotropics are generated in this scattering surface, the existence of 
correlations at angles above about a degree demonstrates that there 
were super-Hubble fluctuations at the recombination epoch. The sign 
of this correlation also confirms the adiabatic paradigm. 

Experimental band-powers for Cf^ from WMAP plus CAPMAP [48], 
CBI [51]. DASI [52], BOOMERANG [53], QUAD [37] and BI- 
CEP [50] are shown in Fig. 23.4. Without the benefit of correlating 
with the temperature anisotropics (i.e., measuring Cj ), the po- 
larization anisotropics are very weak and challenging to measure. 
Nevertheless, there is a highly significant overall detection which is 
consistent with expectation. The QUAD data convincingly show the 
peak at <? ~ 400 (corresponding to the first trough in Cj ) and the 
generally oscillatory structure, while the new BICEP data show the 
lower peak at ^ ~ 140. 

Several experiments have reported upper limits on Cf , but they 
are currently not very constraining. This situation should change as 
increasingly ambitious experiments report results. 

The most distinctive result from the polarization measurements 
is at the largest angular scales (i < 10) in Cj , where there is an 
excess signal compared to that expected from the temperature power 
spectrum alone. This is precisely the signal anticipated from an early 
period of reionization, arising from Doppler shifts during the partial 
scattering at z < Z[. The effect is also confirmed in the WMAP Cf 
results Sit C. = 2-6. The amplitude of the signal indicates that the first 
stars, presumably the source of the ionizing radiation, formed around 
z :^ 10 (somewhat lower than the value suggested by the first year 
WMAP results, although the uncertainty is still quite large). Since 
this corresponds to scattering optical depth r ~ 0.1, then roughly 10% 
of CMB photons were rescattered at the reionization epoch, with the 
other 90% last scattering at 2 ~ 1100. 
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Figure 23.3: Cross power spectrum of the temperature 
anisotropies and E-mode polarization signal from WMAP [47], 
together with estimates from BOOMERANG, DASI, QUAD, 
CBI and BICEP, several of which extend to higher i. Note that 
the widths of the bands have been suppressed for clarity, but 
that in some cases they are almost as wide as the features in tlie 
power spectrum. Also note that the j/-axis here is not multiplied 
by the additional i, which helps to show both the large and small 
angular scale features. 
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Figure 23.4: Power spectrum of E-mode polarization from 
several different experiments, plotted along with a theoretical 
model which fits WMAP plus other CMB data. 

23.7. Complications 

Tliere are a number of issues wliich complicate the interpretation 
of CMB anisotropy data (and are considered to be signal by many 
astrophysicists), some of which we sketch out below. 

23.7.1. Foregrounds : 

The microwave sky contains significant emission from our Galaxy 
and from extra-galactic sources [55]. Fortunately, the frequency 
dependence of these various sources is in general substantially different 
from that of the CMB anisotropy signals. The combination of Galactic 
synchrotron, bremsstrahlung, and dust emission reaches a minimum 
at a wavelength of roughly 3mm (or about 100 GHz). As one moves 
to greater angular resolution, the minimum moves to slightly higher 
frequencies, but becomes more sensitive to unresolved (point-like) 
sources. 

At frequencies around 100 GHz, and for portions of the sky away 
from the Galactic Plane, the foregrounds are typically 1 to 10% of the 
CMB anisotropies. By making observations at multiple frequencies, 
it is relatively straightforward to separate the various components 
and determine tlie CMB signal to the few per cent level. For greater 



sensitivity, it is necessary to use the spatial information and statistical 
properties of the foregrounds to separate them from the CMB. 

The foregrounds for CMB polarization are expected to follow 
a similar pattern, but are less well studied, and are intrinsically 
more complicated. The three year WMAP data have shown that 
the polarized foregrounds dominate at large angular scales, and that 
they must be well characterized in order to be discriminated [56]. 
Whether it is possible to achieve sufficient separation to detect B-mode 
CMB polarization is still an open question. However, for the time 
being, foreground contamination is not a fundamental limit for CMB 
experiments. 

23.7.2. Secondary Anisotropies : 

Witli increasingly precise measurements of the primary anisotropies, 
there is growing theoretical and experimental interest in 'secondary 
anisotropics,' pushing experiments to higher angular resolution and 
sensitivity. These secondary effects arise from the processing of 
the CMB due to ionization history and the evolution of structure, 
including gravitational lensing and patchy reionization effects [57]. 
Additional information can thus be extracted about the Universe at 
z <g; 1000. This tends to be most effectively done through correlating 
CMB maps with other cosmological probes of structure. Secondary 
signals are also typically non-Gaussian, unlike the primary CMB 
anisotropies. 

23.7.3. Sunyaev-Zel'dovich Effect : 

A secondary signal of great current interest is the Sunyaev- 
Zel'dovich (SZ) effect [58], which is Compton scattering (76 — )■ 7'e') 
of the CMB photons by hot electron gas, which creates spectral 
distortions by transferring energy from the electrons to the photons. 
It is particularly important for clusters of galaxies, through which one 
observes a partially Comptonized spectrum, resulting in a decrement 
at radio wavelengths and an increment in the submillimeter. 

The imprint on the CMB sky is of the form AT/T = y f{x), with 
the y-parameter being the integral of Thomson optical depth times 
kT^/nieCr through the cluster, and f{x) describing the frequency 
dependence. This is simply x coth(j:/2) — 4 for a non-relativistic gas 
(the electron temperature in a cluster is typically a few keV), where 
the dimensionless frequency x = hv/kTj. As well as this 'thermal' SZ 
effect, there is also a smaller 'kinetic' effect due to the bulk motion 
of the cluster gas, this being AT/T ~ t{v/c), with either sign, but 
having the same spectrum as the primary CMB anisotropies. 

A significant advantage in finding galaxy clusters this way is that 
the SZ effect is largely independent of redshift, so in principle clusters 
can be found to arbitrarily large distances. The SZ effect can be used 
to find and study individual clusters, and to obtain estimates of the 
Hubble constant. There is also the potential to constrain the equation 
of state of the Dark Energy through counts of detected clusters as a 
function of redshift [59]. Many experiments (including the Planck 
satellite) are currently in operation which will probe clusters in this 
way. The promise of the method has been realised through the first 
detections of clusters purely through the SZ effect, by the SPT [60] 
and ACT [61] experiments. 

23.7.4. Higher-order Statistics : 

Although most of the CMB anisotropy information is contained 
in the power spectra, there will also be weak signals present in 
higher-order statistics. These statistics can measure any primordial 
non-Gaussianity in the perturbations, as well as non-linear growth 
of the fluctuations on small scales and other secondary effects (plus 
residual foreground contamination of course). Although there are an 
infinite variety of ways in which the CMB could be non-Gaussian, 
there is a generic form to consider for the initial conditions, 
where a quadratic contribution to tlie curvature perturbations is 
parameterized through a dimensionless number /nl- This weakly 
non- linear component can be constrained through measurements of the 
bispectrum or Minkowski functionals, for example. Tlie constraints 
depend on the shape of the triangles in harmonic space, and it has 
become common to distinguish the 'local' or 'squeezed' configuration 
(in which one side is much smaller than the other two) from the 
'equilateral' configuration. The results from the WMAP team are 
^9 < /nl < 111 (95% confidence region), for the local mode and 
— 151 < /nl < 253 for the equilateral mode [14]. Different estimators 
used by otlier autliors give results of similar magnitude [62] . 
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The level of /nl expected is small, so that a detection of /j^l ~ 10 
would rule out all single field, slow-roll inflationary models. However, 
with the capabilities of Planck and other future experiments, it seems 
that a measurement of primordial non-Gaussianity may be feasible for 
a wide class of models, and therefore much effort is expected to be 
devoted to predictions and measurements in the coming years. 

23.8. Constraints on Cosmologies 

The most striking outcome of the newer experimental results is that 
the standard cosmological paradigm is in very good shape. A large 
amount of high precision data on the power spectrum is adequately fit 
with fewer than 10 free parameters. The framework is that of FRW 
models, which have nearly flat geometry, containing Dark Matter and 
Dark Energy, and with adiabatic perturbations having close to scale 
invariant initial conditions. 

Within this framework, bounds can be placed on the values of the 
cosmological parameters. Of course, much more stringent constraints 
can be placed on models which cover a restricted parameter space, 
e.g., assuming that fitot = l,n = lorr = 0. More generally, the 
constraints depend upon the adopted prior probability distributions, 
even if they are implicit, for example by restricting the parameter 
freedom or their ranges (particularly where likelihoods peak near the 
boundaries), or by using different choices of other data in combination 
with the CMB. When the data become even more precise, these 
considerations will be less important, but for now we caution that 
restrictions on model space and choice of priors need to be kept in 
mind when adopting specific parameter values and uncertainties. 

There are some combinations of parameters that fit the CMB 
anisotropics almost equivalently. For example, there is a nearly exact 
geometric degeneracy, where any combination of fim and S^a tha-t gives 
the same angular diameter distance to last scattering will give nearly 
identical Ci&. There are also other less exact degeneracies among the 
parameters. Such degeneracies can be broken when using the CMB 
results in combination with other cosmological data sets. Particularly 
useful are complementary constraints from galaxy clustering, the 
abundance of galaxy clusters, baryon acoustic oscillations, weak 
gravitational lensing measurements. Type la supernova distances, 
and the distribution of Lyman a forest clouds. For an overview of 
some of these other cosmological constraints, see The Cosmological 
Parameters — Sec. 21 of this Review. 

The 5-year WMAP data alone, together with constraints from 
Hubble constant determination [63], supernovae [64] and baryon 
acoustic oscillations [29], within the context of a 6 parameter 
family of models (which fixes Htot = 1 and r = 0), yield the 
following results [14]: A = (2.44 ± 0.10) x lO"'', n = 0.960 ± 0.013, 
/i = 0.705 ±0.013, il^h"^ = 0.0227 ± 0.0006, Om/i^ = 0.136 ±0.004 and 
T = 0.084 ± 0.016. There has been little substantive change compared 
with the 3-year data, although the error bars are reduced with more 
integration time and the inclusion of lower frequency polarization data. 
The better measurement of the third acoustic peak, together with 
improved understanding of calibration issues has led to tighter error 
bars on dark matter density and overall normalization. The evidence 
for non-zero reionization optical depth is now very compelling, while 
the evidence for n < 1 is still only at the roughly 3cr level. 

Other combinations of data, e.g., including additional CMB 
measurements, or using cosmological constraints, lead to consistent 
results to those given above, sometimes with smaller error bars, 
and with the precise values depending on data selection [30,37,65]. 
Note that for h, the CMB data alone provide only a very weak 
constraint, unless spatial fiatness or some other cosmological data are 
used. For Q^^h? , the precise value depends sensitively on how much 
freedom is allowed in the shape of the primordial power spectrum 
(see Big-Bang nucleosynthesis — Sec. 20 of this Review). The addition 
of other data-sets also allows for constraints to be placed on further 
parameters. 

For Otot, perhaps the best WMAP constraint is 1.006 ± 0.006, 
from the combination with supernova and baryon acoustic oscillation 
constraints (and setting w = —1). The 95% confidence upper limit 
on r is 0.43 using WMAP alone, tightening to r < 0.22 with the 
addition of other data [14]. This limit depends on how the slope n 
is restricted and whether dn/dlnk ^ is allowed. Nevertheless, it is 
clear that X<l> (sometimes called self-coupled) inflation is disfavored 



by the data, while the m <j> (sometimes called mass term) inflationary 
model is still allowed [14]. Gravity wave constraints coming directly 
from B-mode limits are at the level of r < 0.73 [50]. 

There are also constraints on parameters over and above the basic 
8 that we have described, usually requiring extra cosmological data 
to break degeneracies. For example, the addition of the Dark Energy 
equation of state w adds the partial degeneracy of being able to fit 
a ridge in (w,h) space, extending to low values of both parameters. 
This degeneracy is broken when the CMB is used in combination 
with independent ifg limits, or other data. WMAP plus supernova 
and large-scale structure data yield —0.14 < 1 + w < 0.12 (95% 
confidence), with stronger constraints for flat models. 

For the optical depth r, the best-flt corresponds to a reionization 
redshift centered on 11 in the best-flt cosmology, and assuming 
instantaneous reionization. This redshift appears to be higher that 
that suggested from studies of absorption in high-z quasar spectra [66]. 
The excitement here is that we have direct information from CMB 
polarization which can be combined with other astrophysical 
measurements to understand when the first stars formed and brought 
about the end of the cosmic dark ages. 

23.9. Particle Physics Constraints 

CMB data are beginning to put limits on parameters which are 
directly relevant for particle physics models. For example, there is a 
limit on the neutrino contribution Q^h < 0.0071 (95% confidence) 
from a combination of WMAP and other data [14]. This directly 
implies a limit on neutrino mass, ^ nii, < 0.67 eV, assuming the 
usual number density of fermions which decoupled when they were 
relativistic. Some tighter constraints can be derived using the CMB in 
combination with other data-sets [67]. 

The current suite of data suggest that n < I, with a best-fitting 
value about 5% below unity. If borne out, this would be quite 
constraining for inflationary models. Moreover, this gives a real target 
for B-mode searches, since the value of r in simple models may be in 
the range of detectability, e.g., r ~ 0.2 for m <j) inflation if n ~ 0.95. 
In addition, a combination of the WMAP data with other data-sets 
constrains the running of the spectral index, although at the moment 
there is no evidence for dn/dlnk ^ [14]. 

One other hint of new physics lies in the fact that the quadrupole 
and possibly some of the other low i modes seem anomalously low 
compared with the best-flt ACDM model [6]. Additionally there is 
some weak evidence for a large scale modulation of the smaller-scale 
power [16]. These effects might be expected in a universe which 
has a large-scale cut-off or anisotropy in the initial power spectrum, 
or is topologically non-trivial. However, cosmic variance, possible 
foregrounds, apparent correlations between modes (as mentioned in 
Sec. 23.2), etc., limit the significance of these anomalies. 

In addition, it is also possible to put limits on other pieces 
of physics [68], for example the neutrino chemical potentials, 
contribution of Warm Dark Matter, decaying particles, time variation 
of the fine-structure constant, or physics beyond general relativity. 
Further particle physics constraints will follow as the anisotropy 
measurements increase in precision. 

Careful measurement of the CMB power spectra and non- 
Gaussianity can in principle put constraints on physics at the highest 
energies, including ideas of string theory, extra dimensions, colliding 
branes, etc. At the moment any calculation of predictions appears to 
be far from definitive. However, there is a great deal of activity on 
implications of string theory for the early Universe, and hence a very 
real chance that there might be observational implications for specific 
scenarios. 

23.10. Fundamental Lessons 

More important than the precise values of parameters is what we 
have learned about the general features which describe our observable 
Universe. Beyond the basic hot Big Bang picture, the CMB has 
taught us that: 

• The Universe recombined at z ~ 1100 and started to become 
ionized again at z ~ 10. 

• The geometry of the Universe is close to flat. 



23. Cosmic microwave background 



• Both Dark Matter and Dark Energy are required. 

• Gravitational instability is sufficient to grow all of the observed 
large structures in the Universe. 

• Topological defects were not important for structure formation. 

• There are 'synchronized' super-Hubble modes generated in the 
early Universe. 

• The initial perturbations were adiabatic in nature. 

• The perturbations had close to Gaussian {i.e., maximally random) 
initial conditions. 

It is very tempting to make an analogy between the status of 
the cosmological 'Standard Model' and that of particle physics (see 
earlier Sections of this Review). In cosmology there are about 10 free 
parameters, each of which is becoming well determined, and with a 
great deal of consistency between different measurements. However, 
none of these parameters can be calculated from a fundamental theory, 
and so hints of the bigger picture, 'physics beyond the Standard 
Model,' are being searched for with ever more ambitious experiments. 

Despite this analogy, there are some basic differences. For one 
thing, many of the cosmological parameters change with cosmic epoch, 
and so the measured values are simply the ones determined today, 
and hence they are not 'constants,' like particle masses for example 
(although they are deterministic, so that if one knows their values at 
one epoch, they can be calculated at another). Moreover, the number 
of parameters is not as fixed as it is in the particle physics Standard 
Model; different researchers will not necessarily agree on what the 
free parameters are, and new ones can be added as the quality of 
the data improves. In addition, parameters like r, which come from 
astrophysics, are in principle calculable from known physical processes. 
On top of all this, other parameters might be 'stochastic' in that they 
may be fixed only in our observable patch of the Universe or among 
certain vacuum states in the 'Landscape' [70]. 

In a more general sense, the cosmological 'Standard Model' is 
much further from the underlying 'fundamental theory,' which will 
ultimately provide the values of the parameters from first principles. 
Nevertheless, any genuinely complete 'theory of everything' must 
include an explanation for the values of these cosmological parameters 
as well as the parameters of the Standard Model of particle physics. 

23.11. Future Directions 

Given the significant progress in measuring the CMB sky, which has 
been instrumental in tying down cosmological parameters, what can we 
anticipate for the future? There will be a steady improvement in the 
precision and confidence with which we can determine the appropriate 
cosmological model and its parameters. Ground-based experiments 
operating at smaller angular scales will over the next few years provide 
significantly tighter constraints on the damping tail. New polarization 
experiments will continue to push down the constraints on primordial 
B-modes. The third generation CMB satellite mission, Planck was 
launched successfully in May 2009, and the first results are keenly 
anticipated. 

Despite the increasing improvement in the results, the addition of 
the latest experiments has not significantly changed the established 
cosmological model. It is, therefore, appropriate to ask: what should 
we expect to come from Planck and from other future experiments, 
including those being discussed in the U.S. and in Europe? Planck 
certainly has the advantage of high sensitivity and a full-sky survey. 
A precise measurement of the third acoustic peak provides a good 
determination of the matter density; this can only be done by 
measurements which are accurate relative to the first two peaks 
(which themselves constrain the curvature and the baryon density). A 
detailed measurement of the damping tail region will also significantly 
improve the determination of n and any running of the slope. Planck 
should be capable of measuring Cf quite well, providing both a 
strong check on the cosmological Standard Model and extra constraints 
that will improve parameter estimation. 

A set of cosmological parameters is now known to roughly 10% 
accuracy, and that may seem sufficient for many people. However, 
we should certainly demand more of measurements which describe 
the entire observable Universe! Hence a lot of activity in the coming 
years will continue to focus on determining those parameters with 
increasing precision. This necessarily includes testing for consistency 



among different predictions of the cosmological Standard Model, and 
searching for signals which might require additional physics. 

A second area of focus will be the smaller scale anisotropics 
and 'secondary effects.' There is a great deal of information about 
structure formation at 2 <g 1000 encoded in the CMB sky. This 
may involve higher-order statistics as well as spectral signatures, 
with many new experiments targeting the galaxy cluster SZ effect. 
Such investigations can also provide constraints on the Dark Energy 
equation of state, for example. Planck, as well as new telescopes 
aimed at the highest f.s, should be able to make a lot of progress in 
this arena. 

A third direction is increasingly sensitive searches for specific 
signatures of physics at the highest energies. The most promising of 
these may be the primordial gravitational wave signals in C^ , which 
could be a probe of the ~ 10 GeV energy range. As well as Planck, 
there are several ground- and balloon-based experiments underway 
which are designed to probe the polarization B-modes. Whether the 
amplitude of the effect coming from inflation will be detectable is 
unclear, but the prize makes the effort worthwhile, and the indications 
that n ~ 0.95 give some genuine optimism that r(= T/3) may be of 
order 0.1, and hence within reach soon. 

Anisotropics in the CMB have proven to be the premier probe of 
cosmology and the early Universe. Theoretically the CMB involves 
well-understood physics in the linear regime, and is under very good 
calculational control. A substantial and improving set of observational 
data now exists. Systematics appear to be well understood and not 
a limiting factor. And so for the next few years we can expect an 
increasing amount of cosmological information to be gleaned from 
CMB anisotropics, with the prospect also of some genuine surprises. 
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